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SUMMARY 


In  order  Co  develop  a  model  of  long  nave  atmospheric  noise, 
Pacific-Sierra  Research  Corporation  (PSR)  used  its  full-wave  propaga¬ 
tion  code  to  calculate  waveguide  mode  parameters  in  spread-debris 
nuclear  environments.  Those  parameters  will  be  stored  for  retrieval 
whenever  the  model  is  exercised.  Such  precalculation,  storage,  and 
retrieval  is  far  more  efficient  than  recalculating  mode  parameters 
each  time  the  model  is  used.  A  similar,  but  far  less  extensive  sec  of 
mode  parameters  was  calculated  and  graphed  by  Field  and  Dore  [1975], 
who  treated  only  the  attenuation  rate  of  the  lowest  order  TM  and  TE 
modes  for  propagation  over  sea  water. 

PSR  has,  over  the  years,  received  many  requests  for  those  graphs, 
as  well  as  numerous  inquiries  as  to  (1)  the  behavior  of  parameters 
other  than  attenuation  rate  and  (2)  the  dependence  of  mode  parameters 
on  ground  conductivity.  Because  the  noise-model  data  encompass  the 
parameters  of  all  significant  modes  for  a  wide  range  of  ground  conduc¬ 
tivities,  frequencies,  and  nuclear  environment  intensities,  it  seemed 
worthwhile  to  publish  graphs  of  those  parameters  in  handbook  format. 
For  all  practical  purposes  then,  this  volume  is  hat  handbook. 

The  electron  and  ion  density  profiles  for  the  spread-debris 
nuclear  environments  and  the  collision  frequency  profiles  are  given 
and  discussed  in  Vol.  I  of  the  present  report  and  in  more  detail  bv 
Field  and  Dore  [1975].  The  mode  parameters  plotted  are  attenuation 
rate,  excitation  factor,  phase  velocity,  and  eigencosine.  The  number 
of  modes  considered  varies  from  as  few  as  two,  for  intense  environ¬ 
ments  where  attenuation  is  severe,  to  as  many  as  five,  for  weak  en¬ 


vironments  where  attenuation  is  low  and  more  modes  must  be  retained. 

By  using  the  graphs  in  conjunction  with  Eq .  (6)  in  Vol.  1  and  a  hand¬ 
held  calculator,  estimates  of  nuclear-weapon-induced  signal  loss  can 
be  made . 

This  volume  consists  of  three  sections,  each  of  which  presents 
the  numerical  data  in  a  somewhat  different  form.  Section  1  plots  mode 
parameters  versus  ground  conductivity  for  eight  spread-debris  eiiviron- 
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merits  that  span  the  range  2  x  10"^  >  W  >  2  x  10~^-5  and  frequencies 
that  span  the  range  15  -  45  kHz  in  5  kHz  steps.1'  These  data  make  up 
our  data  base;  Secs.  2  and  3  display  the  data.  Section  2  plots  the 
mode  parameters  versus  W  and  thus  extends  the  graphs  given  by  Field 
and  Dore  [19/5]  .  Section  3  plots  the  mode  parameters  versus  fre¬ 
quency  . 

We  use  the  definitions  of  excitation  values  and  height  gains 
given  by  Pappert  et  al.  [1970]  for  the  Ez  field  from  a  vertical  dipole 
for  transverse  magnetic  (TM)  modes  and  for  the  Ev  field  from  a 
horizontal  dipole  for  transverse  electric  (TE)  modes.  The  TM  height 
gains  are  unity  on  the  ground  and  up  to  aircraft  altitudes;  the  TE 
height  gains  can  be  approximated  by 

fj.(z)  =  n  kz  .  (1) 

5 

Because  TE  mode  attenuation  rates  are  independent  of  ground 
conductivity,  we  do  not  plot  that  data  in  Sec.  1.  However,  data  for 
the  least  attenuated  TE  mode  appears  in  Secs.  2  and  3.  The  TE  ex¬ 
citation  value  is  proportional  to  the  Ey  field  on  the  ground  and  does 
depend  on  the  conductivity.  In  fact,  the  TE  excitation  value  is 
almost  directly  proportional  to  l/n|.  Note  that  the  height-gain 
functions  are  nearly  proportional  to  ng,  so  the  product  of  the  TE 
excitation  value  and  the  two  height-gain  functions  is  virtually  inde¬ 
pendent  of  ground  conductivity  at  aircraft  altitudes. 

The  data  presented  here  were  generated  using  our  anisotropic  full- 
wave  computer  code,  WAVEPROP.  The  WAVEPROP  code  had  its  origins  In  an 
isotropic  code  developed  by  Field,  Lewinstein,  and  Dore  [1976] ,  with 
major  additions  based  on  the  anisotropic  code  developed  at  the  Naval 
Ocean  Systems  Center  (see,  for  example,  Ferguson  and  Snyder,  1987), 
These  codes  use  model  ionospheres  which  specify  the  number  density  of 
electrons  and  ions  at  altitudes  above  the  earth.  They  then  solve  an 


“W  is  the  well-known  spread-debris  Intensity  parameter  [Knapp  and 
Schwartz,  1975;  Field  and  Dore,  1975]. 
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equation,  called  the  modal  equation,  of  the  form 


M(a  ,w,C)  =  0  ,  (2) 

O 

for  the  eigenvalue  C,  which  is  the  cosine  of  the  complex  angle  of 
incidence  of  the  wave  with  the  ionosphere.  In  the  figures  presented 
in  Sec.  1,  the  modes  are  numbered  in  order  of  their  attenuation  rates 
at  high  conductivity  (<7g  =  1  S/m)  .  To  number  the  modes  at  lower 
values  of  <7g,  we  decreased  a g  in  small  steps  and  calculated  mode 
parameters  at  each  step.  Each  mode  was  assigned  the  number  of  the 
mode  that  was  nearest  to  it  (in  C  value)  at  the  prior  step. 

The  fourth  graph  in  each  set  shows  the  path  of  C  in  the  complex  C 
plane.  Section  1  shows  the  path  as  conductivity  changes,  Sec.  2  as  W 
changes,  and  Sec.  3  as  frequency  changes.  As  noted  on  the  graphs,  one 
end  of  the  path  is  marked.  The  graphs  in  Sec.  1  show  both  TE  and  TM 
modes.  Since  C  values  for  the  TE  modes  are  not  sensitive  to  they 
do  not  appear  to  have  a  tail  in  the  plots.  In  certain  cases,  as 
conductivity  decreases,  the  attenuation  rate  for  the  second  TM  mode 
becomes  very  high  and  exceeds  the  attenuation  of  the  third  (and  even 
higher)  modes.  That  behavior  occurs  when  C  approaches  the  cosine  of 
the  Brewster's  angle,  and  we  call  this  mode  the  Brewster's  mode.  The 
Brewster's  mode  is  easily  identified  by  its  path  in  the  c-plane.  In 
Fig.  55,  for  example,  model  2  is  the  Brewster's  mode.  Note  that  its 
path  is  much  longer  than  the  other  modes,  auu  that:  the  path  passes 
over  the  paths  of  several  other  modes.  That  is  when  the  real  part  of 
the  Brewster's  mode  is  the  same  as  the  re.  1  part  of  another  mode,  the 
imaginary  part  of  the  Brewster's  mode  is  larger.  Indeed  the  only  time 
that  two  modes  have  the  same  real  part  is  when  one  of  them  is  a 
Brewster's  mode.  Note  also  in  this  figure  that  the  TM  modes  that  lie 
inside  the  path  of  the  Brewster's  mode  (modes  3  and  4  in  Fig.  55)  loop 
counter  clockwise,  whereas  the  TM  mode  that  lies  outside  it  (mode  5) 
loops  clockwise.  By  looping  counter  clockwise,  we  mean  that  the  real 
part  of  C  decreases  as  conductivity  gets  smaller,  whereas  in.  the  paths 
that  loop  clockwise ,  the  real  part  of  C  increases  as  conductivity 
decreases.  For  the  TM  modes  that  lie  inside  the  Brewster's  mode,  the 
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imaginary  part  of  C  increase  as  conductivity  decreases.  The  maximum 
is  reached  at  the  conductivity  where  the  real  part  of  C  is  the  same  as 
the  real  part  of  the  cosine  of  the  Brewster's  angle. 

Section  1  shows  all  significant  modes;  Secs.  2  and  3  show  only 
the  two  modes  with  the  lowest  attenuation  rate.  In  certain  cases  the 
second  least  attenuated  mode  is  the  second  mode  at  high  conductivities 
and  the  third  mode  at  low  conductivities.  This  mode-switch  causes  a 
discontinuity  in  the  excitation  and  phase  velocity  graphs.  However 
the  program  that  generated  the  figures  uses  a  cubic  spline  to  interpo¬ 
late  between  the  points,  and  so  the  curves  in  those  graphs  appear  to 
oscillate  near  the  transition  point  (see  for  example  Fig.  76).  This 
is  unphysical,  but  the  graphs  were  unreadable  until  the  cubic  spline 
was  used. 

We  define  a  mode  as  being  either  TE  or  TM  by  examining  the  ratio 

„ion  gnd 

1^1  J.K  J.  -  1  /  o  V 

r  =  ~ — r —  •  -  - -  .  (3) 

,,R;on,,ard  - 1 


The  numerator  and  denominator  of  this  equation  are  the  isotropic  TE 
and  TM  modal  equations,  respectively.  Recall  that  the  condition  for  a 
mode  is  that  the  model  equation  be  zero.  Therefore,  r  »  1  for  TM 
modes,  and  r  «  1  for  TE  modes.  For  the  disturbed  ionospheres  con¬ 
sidered  here,  it  was  easy  to  identify  and  label  the  type  of  mode. 

Under  nighttime  conditions,  the  ratio  [Eq.  (3)]  might  approach  one 
which  would  make  a  clear  distinction  difficult. 

The  data  in  this  volume  were  generated  assuming  that  the  electron- 
neutral  collision  frequency  was; 

t>(z)  =  1.8  x  10  e  '  s  .  (4) 

Here  z  is  the  altitude  in  kilometers.  The  mass  of  the  ions  is  taken 
to  be  32  AMU,  and  the  ion-neutral  collision  frequency  is  1/40  of 
Eq.  (4).  The  relative  permittivity  of  the  ground  is  taken  to  be  15. 
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CONVERSION  TABLE 


[ 

f  Conversion  factors  for  U.S.  Costomary  to  metric  (SI)  units  of  measurement 


[  MULTIPLY - - 

BY 

*>  TO  GET 

1  to  git;  <4 - - - 

i 

-  BY  + - 

-  DIVIDE 

Angstrom 

1. 1)00  000  X  E  -10 

meters  (mi 

|  atmosphere  (normal) 

1  C1J  25  X  E  +2 

kilo  pascal  (kPaj 

k  b*r 

1.  000  000  X  E  *2 

kilo  psecal  (kPa) 

i  hum 

1.  000  000  X  E  -23 

•7  •% 

rreter"  rm") 

►  British  thermal  unit  (thermo  Chemical) 

X.  0&4  350  X  £  +3 

joule  (J) 
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4. 184  000 

joule  (J) 

(  cal  (thermochemical) /cm" 

4.  184  000  X  E  -2 

mega  joule/m"  {MJ/m") 

curie 

3 .  TOO  00r  X  E  -  1 

•gigs  bocquarel  (GBq) 

1  degree  ianglfc) 

1.  ?4$  319  X  E  -2 

radian  (rad) 

„  degree  Fahrenheit 

<«  •  (t*  f  *  459,  671/1.  8  • 
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1.  602  13  X  E  *19 

joulu  (3) 

;  erf[ 
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1.  000  000  X  E  -7 
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3.  048  000  X  E  -1 

mete  j  (m) 

1  f  co  t -pound  -force 
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joule  (J) 

^  gallon  [L1.  S.  liquid) 

3.  785  412  X  E  -3 

Mstsr3  (m3) 

inch 

2.  V0  000  X  E  -2 

mot«r  lm) 

[  jork 

1.  000  000  XE  *9 

joule  (J) 

!  jouls/Hilesram  iJ/Tg)  1  radiation  dose 

j  absortwdj 

1.000  000 

Gray  (Gy) 

^  kilotona 

4.  183 

terajoulas 

1  kip  (1000  Ibt) 

4.  448  222  X  E  *3 

naar  ton  (N) 

,  kip/lnch2  (ksi) 

S.  994  757  X  E  *3 

kilo  pascal  (kPa) 

l  kup 

nawton  -Me  end  /m3 

[ 

1.000  000  X  E  *2 

(N-a/m2) 

micron 

1  000  000  X  £  45 

motor  (m) 

mil 

2.  540  000  X  E  -5 

ma'ar  (m) 
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1.  609  344  XE*3 
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4.  448  222 
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pound -force  uich 
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pound-mass  (Ibm  avoirdupois) 

4.  535  924  X  E  -1 

kilogram  (kg) 

pound -mass -foot 2  (moment  of  inertia) 

kilogram -molar3 

4. 214  011  X  E  -2 

Og-m2) 

pound -mu«  /Toot3 

3 

kilogram  /meter 

1.  601  845  X  E  ♦  ! 

(kg/m'-z 

rmd  (radiation  dose  absorbed) 
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••Gray  (Gy) 

roentgen 

emilomb  /kilogram 

: 
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(C/kg) 
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second  (a) 
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torr  (mm  H*.  0*  C) 
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■The  b*cqusrel  (Bq)  n  (ha  SI  unit  of  radioactivity;  1  Bq  .  1  tvant/a. 
’■'Die  Gray  (Gy)  it  the  31  unit  of  absorbed  radiation. 
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Figure  2.  Parameters  for  W  »  2  x  1 0’®,  frequency  -  20  KHz. 


4 


<  -40 

-20 

10'5  10'4  10'3  10‘ 

Cg  (S/m) 


b.  TM  excitation  values  as  function  of  ground  conductivity. 
Mode  1  — —  Mode  2  . 
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Figure  4.  Parameters  for  W  «  2  x  10'®,  frequency  -  30  kHz. 
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a.  Attenuation  as  function  of  ground  xnductivity 
Mode  1  -  Mode  2 . 


Figure  5.  Parameters  for  W  «  2  x  1 0‘®,  frequency  -  35  kHz. 


a.  Attenuation  as  function  of  ground  conductivity 
Mode  1  -  Mode  2  . 
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8,  Parameter?  for  W  -  2  x  10'®,  frequency  »  15  KHz. 
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a.  Attenuation  as  function  of  ground  conductivity 
Mode  1  -  Mode  2  . 
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b.  TM  excitation  values  as  function  of  ground  conductivity. 
Mode  1  -  Mode  2 . - 
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Figure  9.  Parameters  for  W  -  2  x  1 0'9,  frequency  -  20  kHz. 
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a.  Attenuation  as  function  of  ground  conductivity 
Mode  1  -  Mode  2 . 


b.  TM  excitation  values  as  function  of  ground  conductivity 
Mode  1  -  Mode  2 . 
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b.  TM  excitation  values  as  function  of  ground  conductivity. 
Mode  1  - -  Mode  2 . 
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Figure  17.  Parameters  for  W  -  2  x  10’10,  frequency  -  25  kHz. 


34 


a.  Attenuation  as  lunctlon  of  ground  conductivity 
Mode  1  -  Mode  2  . 
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Figure  18.  Parameters  for  W  -  2  x  10’10,  frequency  -  30  kHz. 
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a.  Attenuation  as  function  of  ground  conductivity. 
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b.  TM  excitation  values  as  function  of  ground  conductivity. 
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Figure  32.  Parameters  for  W  -  2  x  10'12,  frequency  -  30  kHz. 


c.  Relative  phase  velocity  as  a  function  of  ground  conductivity 
Mode  1  -  Mode  2  .  Mode  3 - 


d.  Mode  paths  in  the  C-plane  as  conductivity  changes 
Mode  1  -  Mode  2 .  Mode  3 - 


NOTE:  High  oondurtivity  point  mwfcwf  wWi  ■  fw  TM  rood**,  A  for  TE. 

Figure  33  Parameters  forW-2  x  10'12,  frequency  -  35 KHz  (Concluded) 


c.  Relative  phase  velocity  as  a  function  of  ground  conductivity. 

Mode  1  -  Mode  2  .  Mode  3 - Mode  4  — —  Mode  5 


d.  Mode  paths  in  the  C-plane  as  conductivity  changes. 
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NOTE;  High  conductivity  point  marked  with  ■  for  TM  mod**,  A  for  Tfc. 

Figure  35.  Parameters  for  W  -  2  x  10'^  frequency  ~  45  KHz  (Concluded) 


Altenualion  rale  (dB/Mm) 
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b.  TM  excitation  values  as  function  of  ground  conductivity. 
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Figure  41 .  Parameters  for  W  -  2  x  1 0'1 3,  frequency  -  40  kHz. 
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a.  Attenuation  as  function  of  ground  conductivity. 
Mode  1  - -  Mode  2  .  Mode  3 - 


b.  TM  excitation  values  as  function  of  ground  conductivity. 
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Figure  47.  Parameters  for  W  ■  2  x  1 0‘14,  frequency  -  35  KHz. 
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c.  Relative  phase  velocity  as  a  function  of  ground  conductivity 
Mode  1  -  Mode  2  .  Mode  3 - 


d.  Mode  paths  in  the  C-plane  as  conductivity  changes 
Mode  1  — —  Mode  2  .  Mode  3  — — 


NOTE:  High  conductMy  point  m*rf<*d  with  V  for  TM  mod**,  a  lor  TE, 

Figure  48.  Parameters  for  W  ■  2  x  10"14,  frequency  -  40  KHz  (Concluded) 


51 .  Parameters  for  W  -  2  x  1 0'1  frequency  -  20  kHz. 
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SECTION  2 


PARAMETERS  AS  A  FUNCTION  OF  W 
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a.  Attenuation  as  (unction  of  parameter  W. 

15  kHz  -  25  kHz  .  35  kHz - 45  kHz 


Figure  61.  Parameters  for  second  least  attenuated  TM  mode,  Gg  -  10 


Attenuation  rate  (dB/Mm) 


Figure  63,  Parameters  for  second  least  attenuated  TM  mode,  Og  ~  3  x  1 0'4  S/m, 
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45  kHz 


a.  Attenuation  as  function  of  parameter  W. 


Figure  67.  Parameters  for  second  least  attenuated  TM  mode,  Cg  -  3  x  ICr5  S/m, 
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PARAMETERS  AS  A  FUNCTION  OF  FREQUENCY 
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Excitation  values  fd£l)  Attenuation  rate  (dB/Mm) 


a.  Attenuation  as  function  of  frequency. 
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b.  TM  excitation  values  as  function  of  frequency. 


Altenualion  rate  (dB/Mm) 


Excitation  values  (dB)  Attenuation  rate  (dB/Mm) 


v/c  - 


c.  Relative  phase  velocity  as  function  of  frequency. 
W«2x10-8  -  W«2x10*10  .  W-2x1(H2 
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d.  Mode  paths  in  the  C-plane  as  frequency  changes. 
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Figure  73.  Parameters  for  least  attenuated  TM  mode,  Og  -  10*3  S/m  (Concluded) 
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b.  TM  excitation  values  as  function  of  frequency. 
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Figure  75.  Parameters  for  least  attenuated  TM  mode,  <jg  -  3  x  10*4  s/m. 
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Figure  76.  Parameters  for  second  least  attenuated  TM  mode,  <jg  -  3  x  1 0'4  S/m 
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b.  TM  excitation  values  as  function  of  frequency. 
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Figure  78.  Parameters  for  second  least  attenuated  TM  mode,  Og  -  10‘4  S/m, 
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c.  Relative  phase  velocity  as  function  of  frequency. 
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d.  Mode  paths  in  the  C-plane  as  frequency  changes. 
W«2x10-8 -  W-2x10-10  .  W-2X10'12  — 


—  W-2x10‘15 - 
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Figure  78.  Parameters  for  second  least  attenuated  TM  mode,  Og  -  10"4  S/m  (Concluded) 
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b.  TM  excitation  values  as  function  of  frequency. 

W-2x10'8  -  W-2x10-10  . .  W-2X10'12 - W-2x1CH5 - 


Frequency (kHz) 

Figure  80.  Parameters  for  second  least  attenuated  TM  mode,  <jg  -  3  x  10"5  S/m. 
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c.  Relative  phase  velocity  as  function  of  frequency. 
W-2x10-8  -  W-2xlO-10  .  W«2x1CH2 


W«2x1(H5 


f-requency  (kHz) 

d.  Mode  paths  in  the  C-piane  as  frequency  changes. 
W-2x10-8  -  W-axICHO  . .  W«2x10'12  - 

-  W-2x10-15 

NOTE;  Thfi  point  for  kwwrt  frequency  maiktd  with  ■ 

Figure  81 .  Parameters  for  least  attenuated  TM  mode,  Og  -  10'5  S/m  (Concluded) 
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